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Abstract. We derive a new criterion for a real-valued function u to be in the 
Sobolev space W 1 ' 2 (R n ). This criterion consists of comparing the value of a 
functional f f(u) with the values of the same functional applied to convolutions 
of u with a Dirac sequence. The difference of these values converges to zero 
as the convolutions approach u, and we prove that the rate of convergence 
to zero is connected to regularity: u £ W 1,2 if and only if the convergence 
is sufficiently fast. We finally apply our criterium to a minimization problem 
with constraints, where regularity of minimizers cannot be deduced from the 
Euler-Lagrange equation. 



1. Introduction 

Jensen's inequality states that if / : K — ► R is convex and <p G L 1 (K") with 
tp > and j tp = 1, then 



f(u(x))(p(x) dx > / I / u(x)ip(x) dx 



for any u : M. n — ► K for which the integrals make sense. A consequence of this 
inequality is that 

f(u) = f f{u)*y > f f(u*<p). (1.1) 

In this paper we investigate the inequality (jl.ip more closely. In particular we 
study the relationship between the regularity of a function u E L 2 (W l ) and the 
asymptotic behaviour as e — > of 

T]{u) := / [/(«)- /(u*v> e )]. (1.2) 

Here <p s {y) — £~ n tp{£~ 1 y) an d / is a smooth function, but now not necessarily 
convex. Since u * (p e — > u almost everywhere, we find for 'well-behaved' / that 
Tf(u) — > as e — > 0. Our aim is to establish a connection between the rate of 
convergence of Tf(u) to zero and the regularity of u. 
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Such a connection between the decay rate of Tj{u) and the regularity of u is 
suggested by the following informal arguments. Taking n = 1, and assuming tp to 
be even and u to be smooth, we develop u * ip E as 



with c = i J R y 2 <p(y) dy. This suggests that if it 6 W 1,2 (W) then Tj (it) is of order e 2 . 
Moreover, since /" > in the case that / is convex, (|1.4[) gives an enhanced version 
of (Ti~Tjl . 

Conversely, for functions not in W 1,2 (R ra ) we might not observe decay of T^(u) 
with the rate of e 2 , as a simple example shows. Take f(u) = u 2 and consider a 
function with a jump singularity such as u(x) — H{x) — H(x — 1) ^ W Al ' 2 (M), where 
H (x) is the Heaviside function. Now choose as regularization kernels the functions 
(fi £ (x) — ^(H(x + e) — H(x — e)), after which an explicit calculation shows that 
T|(m) = 2e/3. Here the decay is only of order e. 

The goal of this paper is to prove the asymptotic development (|1.4p in arbitrary 
dimension and to show that also the converse statement holds true: if T^(u) is 
of order e 2 then u € W 1,2 (W l ). These results are stated below and proved in 
Sections [2] and [3] 

The fact that one can deduce regularity from the decay rate is the original 
motivation of this work. In Section 2] we illustrate the use of this result with a 
minimization problem in which the Euler-Lagrange equation provides no regular- 
ity for a minimizer u. Instead we estimate the decay of T^(u) directly from the 
minimization property and obtain that u is in W ' . 

Our results can be compared to the characterisation of Sobolev spaces introduced 
by Bourgain, Brezis and Mironescu [2j [4]. In fact, the regularity conclusion in 
Theorem 11.21 could be derived from 2J , with about the same amount of effort as 
the self-contained proof that we give here. 

1.1. Notation and assumptions. Let ip £ L 1 (R") satisfy 




(1.3) 



so that 




(1.4) 



ip > in 1", 



(1.5) 




(1.6) 



(1.7) 



(1.8) 



For s > we define the Dirac sequence 




(1.9) 
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Eventually we will restrict ourselves to the case that ip is rotationally symmetric, 
that is <p{x) — (p{\x\), where (p : M.q — > Mq~. Then (11. 6| is equivalent to 



nu n / r"-VeW^ = 1, (1-10) 
Jo 

where lu u denotes the volume of the unit-ball, i.e. |i?i(0)| = u) n . 

For a function u G L x (R n ) and < s < 1 we define the convolution u £ and 
modified convolutions u £ s by 

„.(.)-(..«.)«-/' «.-,*.<»>*, 



u e , s (x) := (u*95 es )(x) = / u(x - y)tp £S (y)dy = / u(x - sz)tp e (z) dz. 

JR™ JR™ 

(1.12) 

Note that u £ = u £ .i,u — it e .o- 

We also use the notation a <g> b for the tensor product of a, b G M. n . 



1.2. Statement of main results. Our first result proves (|1.4[) . 

Theorem 1.1. Lei / G C 2 (K) /iaue uniformly bounded second derivative and 

/(0)=0, /'(0) = 0. (1.13) 
Lei (iy9 e ) e >o be a Dirac sequence as in Section \l.l\ Then for any u G W ' (M™), 

lim-W [f(u)-f(u s )]dx=~[ f"(u(x))Vu(x)-A(<p)Vu(x)dx, (1.14) 

w/iere 

A(y>) = / (y®y)p{y)dy. (1.15) 

JR" 

//<,£> is rotationally symmetric, i.e. p{x) = <p(\x\), then 

1™, p _£ n [/(«) - /(«*)] = J f"(u(x))\Vu(x)\ 2 dx, (1.16) 

a(<p) = ui n / r n+1 cp(r)dr. (1.17) 
Jo 

The second theorem shows that for uniformly convex / a decay of Tf(u) of 
order e 2 implies that u G H /1 ' 2 (R"). 

Theorem 1.2. Let f G C 2 (M) have uniformly bounded second derivative, assume 
that (|1.13[) is satisfied, and that there is a positive number c\ > such that 

f">ci onR. (1.18) 

If 'for ue L 2 {W l ), 

liminfi / [/(u(x)) - /(u e (x))] dx < oo, (1.19) 

£ Jr" 

then u G W 1,2 (R n ). In particular (fTTT4|) and (TTTd]) ftoZd. 

Remark 1.3. We prescribe f|l . 13[) since in general the difference f(u) — f{u £ ) need 
not have sufficient decay to be Lebesgue integrable. For functions u G L 1 (R") R 
L 2 (R n ) Theorems O and O hold even without assuming ([TTT3]l . 
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2. Proof of Theorem 11.11 



We first show that f(u) - f(u e ) 6 L 1 (R n ). Let i) e C°(R"; [0,1])- Using the 
Fundamental Theorem of Calculus we obtain that 



r](x)\f(u(x)) - f(u £ (x))\ dx = 

dx 



f 1 d 

V(x) J ^/(u £)S (x)) ds 



< 



l 

dy da; ds 



^O) / f'{u e ^{x))Vu{x-sy)-yLp E (y) 



-2i I v<yX ^ / ^'( u ^- s ^ x )) 2 \y\^^y) d y dxds 



T)(x) \Vu(x - sy)\ 2 \y\ip £ (y)dydxds, (2.1) 



1 

+ 2 

where we have used Fubini's Theorem and Young's inequality. For the first term 
on the right-hand side we deduce by (|1.8p and (|1.13|) that 



r,(x)f'(u e , s (x)Y\y\ip e (y)dydx < C e (^)[[/"||L / u E<s {xf dx 



< C e { 9 )\\f\\lo J ^{xfdx, ( 2 - 2 ) 

For the second term on the right-hand side of (|2.ip we obtain similarly 

??(a;)|Vu(a;-sy)| 2 |y|^ e (y)dj/da; < C e (</>) / \Vu{x)\ 2 dx. (2.3) 

By ([2~T j) -(|2T3" l) we therefore obtain 

V {x)\f{u{x))-f{u £ {x))\dx < C e (<p)(l + \\f"\\l)\\u\\ 2 wl , 2{Mn) < oo. (2.4) 

Letting 77 /" 1 we deduce that f(u) — f(u e ) G i 1 (M"). Repeating some of the 
calculations above, we obtain that 

TUu) = - f [ ^-f(u £ , s (x))dsdx 

JR" JO OS 

mf {u £ . s (x))Vu(x - sy) ■ yip £ {y)dydxds (2.5) 
_:™ 

and that 

f'(u e , s (x))-f'(u e , s (x-sy)) = - I — f '{u Sj8 {x - ry)) dr 



dr' 

f"(ue }S (x - ry))Vu £ , s (x - ry) -ydr. (2.6) 
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Since for all < s < 1 

f'(u e . s (x - sy))Vu(x - sy) ■ yip s (y) dy dx 

I / /'(«s,iW)Vu(z) • y(p e (y)dzdy 

we deduce from ([23]) . (|2~fif that 

[/(u(i)) - /(u e (i))] dx = 



( / f"(u E , s {x - ry))Vu s<s (x - ry) -ydrj 



T n+2 



"1 />s 
'0 ./0 



Vtt(x — sy) • y<p £ (y) dydxds 

f" (u SjS (x))Vu S!S (x) ■ y 
Vu(i — ry) ■ yip £ (y) dy dx dr ds 
f"(u StS (x))Vu e , s (x)- 

[y ® y) ' Vw(x — ry) (/3(e _1 y) cfa; dr ds 



1 



f" (u e>s (x))Wu EiS (x) 
1 



r n+2 



(z z) ■ Vm(x — z) ip(e r 1 z)dzdxdrds 



(2.7) 



(2.8) 



f" (u s ,s(x))Vu e ,s[x) ■ («w * Vu)(as) dxdr ds, 
with the modified convolution kernel n Er defined by 

Ker(z) := (er) - "*/— V 
where k is given by 

«(y) := {y<E>y)(f(y). (2.9) 
As e — > we have that for all < r, s < 1, 

u e , s — > it almost everywhere, 

in L 2 {W l ), (2.10) 

K(y)dy}Vu inL 2 (R n ). (2.11) 



Vu e , s — > Vu 

Moreover, the integrand on the right-hand side of (|2.7|) is dominated by the function 



^II/"IIl-(e)(|Vu E:S | 2 + \ Ker * Vi*| 2 ), 



(2.12) 
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which converges strongly in L 1 (M™) as e — > by ()2.10j) and (|2.11|) . We therefore 
deduce from (12.711 that 



1 



^2 /_ [f( U ( X )) ~ f( u e(x))] dx 



7* 

Jo 



f"(u(x))Vu(x) ■ ( / K,(y) dy)V u(x) dx dr ds 



= \ j n f"(u(x))Vu(x) • ( j n(y) dyj Vu(x) dx. 
In the rotationally symmetric case we observe that 

(y ®y)ip{y)dy = ( j r n+1 (p(r) dr^j j 6 ® 6 d0 



(2.13) 



(2.14) 



Here we used that for i, j = 1, . . . , n, 



OiQj d6 



u n if i=j, 
if z ^ J. 



We therefore obtain from (|2.7p ~ (|2.14p that in the rotationally symmetric case 



lim 1 J [f(u) - f(u E )] dx = uj n (J~ ' r n+1 f>(r)dr^J^f"(u)\Vu\ 2 dx. 



(2.15) 



3. Proof of Theorem 11.2 



Let 



A := liminf — 

e^O e 2 



[f(u(x)) - f(u £ (x))] dx. 



We first remark that by assumption (JXTTHJ) the function r h /(r) 
and we deduce again by (jl.ip that 



ci 2 
2 £ 



> 0. 



[(/(«) -y- 2 )-(/(«e) 
This implies that 

liminf \ [ (u 2 -uf) < - liminf \ j \f(u)-f(u e )] 



(3.1) 



^-r 2 is convex 



1a. 

ci 



(3.2) 



(3.3) 
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Next we consider S > and the regularizations us — u * (p$. Set 7 £ := (p £ * tp E 
Then we obtain that 

[ug(x) - (u s * <p e ) (x)] dx = 

.2 



u s - / / u s (x)u 5 (y)~f e (x ~ y)dydx 

Iff 2 

(u s (x) -u s (y)) 7 e (x -y)dydx 



2 

ill. Ill \ f f f 2 

- 2J J J ( u ( x - z )- u (y~ z )) fs(z)dzj E (x-y)dydx 
= \J J J M ( u (0- u {v)) 2 Je{£-v)dr]d£,<ps{z)dz 
= \\ I (u(0-u(r,)) 2 lE (li-r,)dr,d£ 

L JR" JR« 

K(0-« 2 (oR- (3.4) 

We therefore deduce from Theorem 11.11 and (|3.3|) , (|3.4|) that 

/ \Vu s \ 2 = C(y)liminf \ [ [u 2 s (x) - (u 5 * ip e ) 2 (x)] dx 

< C(^)liminfi / [w 2 (^) - u 2 e (x)] dx < —AC{tp). (3.5) 

£ Jr" ci 

Since this estimate is uniform in 6 > it follows that u £ W 1 ' 2 (M. n ). By Theo- 
rem [TTT] we therefore deduce f| 1 . 14[) and [)1.16p . 

4. Application: regularity of minimizers in a lipid bilayer model 



We illustrate the utility of Theorems 11.11 and 1 1 . 21 with an example. The problem 
is to determine the regularity of solutions of the following minimization problem. 

Problem 4.1. Let a > 0, h > and a kernel n £ VK 1 ' 1 (R, Kg") be given such that 
k' £ BV(M). Denote by Th the translation operator, 

{T h u)(x) := u(x - h) for w : K ^ R. 

Consider the set K C K, 

K := iu 6 i x (K) I u = l,ti>0,tt|T A tt<lo.e.) (4.1) 
Jr ' 

and a strictly convex, increasing, and smooth function f with /(0) = /'(0) = 0. 
We then define a functional F : K — > K by 

F{u) := I f{u) - I aun*u (4.2) 
Jr. Jr 

and consider the problem of finding a function u £ K that minimizes F in K , i.e. 

F(u) = min{F(u) v £ K}. (4.3) 
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Remark 4.2. This problem arises in the modelling of lipid bilayers, biological 
membranes (see p] for more details). There, a specific type of molecules is consid- 
ered, consisting of two beads connected by a rigid rod. The beads have a certain 
volume, but the rod occupies no space. 

The beads are assumed to be of sub-continuum size, but the rod length is non- 
negligeable at the continuum scale. In the one-dimensional case we also assume 
that the rods lie parallel to the single spatial axis. Combining these assumptions 
we model the distribution of such molecules over the real line by a variable u that 
represents the volume fraction occupied by leftmost beads. The volume fraction of 
rightmost beads is then given by t^m, and in the condition u + t^u < 1 we now 
recognize a volume constraint. 

The functional F represents a free energy. In the entropy term J f(u) the func- 
tion / is strictly convex, increasing, and smooth and can be assumed to satisfy 
/(0) = /'(0) = 0. The destabilizing term — J uk*u is a highly stylized representa- 
tion of the hydrophobic effect, which favours clustering of beads. 

Without the constraints u > and u + t^u < 1 present in (|4.1[) . we would 
immediately be able to infer that minimizers (even stationary points) are smooth 
using a simple bootstrap argument: since u G K we have u G L P (M.) and k*u G W 1,p 
for all 1 < p < oo, and therefore, using the Euler-Lagrange equation, 



for some A G R, we find u G W 1:P (R). Iterating this procedure we obtain that u G 
W k ' p (R) for all k G N , 1 < p < oo. However, when including the two constraints, 
two additional Lagrange multipliers (i and v appear in the Euler-Lagrange equation, 



Here fj, and v are measures on R, and standard theory provides no further regularity 
than this. In this case the lack of regularity in the right-hand side interferes with the 
bootstrap process, and this equation therefore does not give rise to any additional 
regularity. 

The interest of Theorem 11.21 for this case lies in the fact that K is closed under 
convolution. In particular, if u is a minimizer of F in K, then u e — tp £ *u is also 
admissible, and we can compare F(u e ) with F(u). From this comparison and an 
application of Theorem 1 1.21 we deduce the regularity of u: 

Corollary 4.3. Let u be a minimizer of (I4.3[) . Then u G Vt /1:2 (R). 

Proof. Choose a function ip G i 1 (R) as in Section [TTT1 with Dirac sequence (ip e ) £> o, 
and set u e := ip e *u. First we show that there exists a constant C G R such that for 



f'{u) — 2an*u = A, 



/ ' (ll) — 2otK*U = A + /i — V — T-yJJ. 



all u G L 2 (R) 




(4.4) 



With this aim we observe that 




(4.5) 



with 

7 := le{x) := (<p e *tp e )(x) = -7 

e 

and that 7,7 e satisfy the assumptions in Section ITTTI 
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Since k! G BV(R) we have (k * u)" = k" * u G L 2 (M) and 



l K "l)ll w IU 2 (R)- 



Repeating some arguments of Theorem 11.11 we calculate that 

r 1 d 



ds 



(k*u)(x - sy)^ E {y)dyds 



(k*u)'(x - sy)y-y e (y) dy ds 

(k*u)" (x — rsy)sy 2 j E (y) dy ds dr 
(k*u)" ' (x — ersz)sz 2 j(z) dz dr ds, 



and therefore that 



< 



(k*u)" (x — ersz) 2 dx / u(x) 2 dx sz 2, -f(z) dz dr ds 



< £ 2 C7( 7 )|| K "* u || i2(R) || w || 
<e 2 C((p,K)\\u\\ 2 L2(R) , 



which by (|4.5p implies (|4.4|) . 

Since K is closed under convolution with ip e , u e is admissible, and therefore 

< F(u e ) - F(u) < f {f(u e ) - /(«)] + Ce 2 f u 2 . 

The last term is bounded by 

Ce 2 IMI L i (R) \H Laom <Ce 2 , 

which follows from combining both inequality constraints in (|4.ip . Therefore 
f[f(u) - f(u £ )} = 0{e 2 ), and from Theorem O we conclude that u G W 1,2 (R). 

□ 

As is described in more detail in the thesis [3] , Corollory 14.31 paves the way 
towards rigorously deriving the Euler-Lagrange equations for Problem (|4.3|) . and 
forms an important ingredient of the proof of existence of minimizers. 
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